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Evolution of non-linear cosmological perturbations 

David Langlois^’^, Filippo Vernizzi^ 

^APC (Astroparticules et Cosmologie), UMR 7164 (CNRS, Universite Paris 7, CEA, 
Observatoire de Paris), 11 Place Marcelin Berthelot, F-75005 Paris, France 
^Institut d’Astrophysique de Paris, 98bis Boulevard Arago, 75014 Paris, France, and 
^Plelsinki Institute of Physics, P.O. Box 64, FIN-00014, University of Helsinki, Finland 

We define fully non-perturbative generalizations of the uniform density and comoving curvature 
perturbations, which are known, in the linear theory, to be conserved on sufficiently large scales 
for adiabatic perturbations. Our non-linear generalizations are defined geometrically, independently 
of any coordinate system. We give the equations governing their evolution on all scales. Also, in 
order to make contact with previous works on first and second order perturbations, we introduce a 
coordinate system and show that previous results can be recovered, on large scales, in a remarkably 
simple way, after restricting our definitions to first and second orders in a perturbative expansion. 


The relativistic theory of cosmological perturbations is 
a cornerstone in our understanding of the early universe 
as it is indispensable in relating the scenarios of the early 
universe, such as inflation, to cosmological data such as 
the Cosmic Microwave Background (CMB) anisotropies. 
This is the reason why the theory of linear cosmologi¬ 
cal perturbations has been developed to a high de^ee of 
sophistication during the last twenty-five years 0, 0 ■ 

With the increase of precision of the CMB data, the 
study of relativistic cosmological perturbations beyond 
linear order is becoming a topical subject, especially to 
study primordial non-Gaussianity. This is why much ef¬ 
fort has been devoted recently to the investigation of 
second-order perturbations, mainly resorting to an ex¬ 
pansion, up to second order, of Einstein’s equations ii 
This formalism is rather difficult, due to the complexity 
of Einstein’s equations, and the identification of gauge 
invariant quantities is not straightforward 00 ■ 

Other recent approaches of non-linear perturbations 
i II E0 are based on the long wavelength approx¬ 
imation, which at lowest order is related to the so 
called separate universe picture that represents our uni¬ 
verse, on scales larger than the Hubble radius, as juxta¬ 
posed Friedmann-Lemaitre-Robertson-Walker (FLRW) 
universes with slightly different scale factors. 

In the present work, we propose a different approach 
based on a purely geometrical description of the pertur¬ 
bations. Our approach does not rely on any perturbative 
expansion and is thus not limited to second-order per¬ 
turbations. It does not require any approximation and 
covers both large and small scales. Furthermore, as we 
shall see, it can be related rather easily to the formula¬ 
tions mentioned above. Our approach is directly inspired 
by the so called covariant formalism for cosmological per¬ 
turbations introduced by Ellis and Bruni El- At linear 
order, the covariant formalism is computationally equiv¬ 
alent to the much more used coordinate approach. How¬ 
ever, as we show here, it is tremendously more powerful 
when one goes beyond linear order. 

In particular, we define geometrically the generaliza¬ 
tions of quantities widely used in the linear theory be¬ 


cause they are conserved on large scales. We then give 
the full non-linear equations that govern these quantities. 
To emphasize the efficiency of our approach with respect 
to the more traditional coordinate-based approach, we 
show how second order gauge invariant quantities, re¬ 
cently derived in the literature, can be recovered, on large 
scales, from our geometric quantities via a straightfor¬ 
ward derivation. 

Covariant approach - We consider a spacetime, with 
metric gab, filled with a single perfect fluid characterized 
by a four-velocity = dx°‘/dr {uaU°' = —I), where r 

is the proper time along the flow lines, proper energy 
density p and pressure P. The corresponding energy- 
momentum tensor is T‘\ = [pP) u’^Ub + Pg\- We also 
introduce the volume expansion 0 = \IaU°‘, which re¬ 
duces to the Hubble parameter iJ, 0 = iH , in a FLRW 
spacetime. The integration of 0 along the fluid world 
lines, with respect to the proper time r. 



(0 = 3d = 3w“Vaa), 


( 1 ) 


can be used to define, for each observer comoving with 
the fluid, a local scale factor e“. Note that a is defined 
up to an integration constant for each fluid world line, 
which we can set to zero on some reference hypersurface. 

As shown by Ellis and Bruni 0, it is possible to de¬ 
fine, in a geometrical way, quantities that can be inter¬ 
preted as perturbations with respect to the FLRW con¬ 
figuration, by introducing spatially projected gradients. 
The spatial projection tensor, orthogonal to the fluid ve¬ 
locity m “, is defined by h\ = g\ -\- u’^Ub- Defining the 
spatial gradient Da = haVb, we consider the quantities 


Xa = DaP, Ya = DaP, Za = DaQ, Wa = Da 


which automatically vanish in a strictly FLRW space- 
time: in this sense we call them perturbations. However, 
they are fully non-perturbative quantities. Note that the 
last quantity, Wa, which as far as we know has not been 
introduced previously in the covariant formalism, plays a 
crucial role in our approach, where we replace the famil¬ 
iar spatial curvature perturbation by the perturbation of 


2 


the integrated expansion a, i.e., of the local number of 
e-fo lds, in accordance with the separate universe picture 

dlla 

Evolution of perturbations - Starting from the above 
definitions, we now derive the fundamental evolution 
equation for these non-linear perturbations. As in m 
in the context of the linear theory, our starting point 
here will be the energy-momentum conservation, 

v^n = 0. (2) 

The projection along of Equation 0 yields 
p + Q{p + P) = 0, 

where a dot represents a derivative with respect to r, i.e., 
stands for u°‘Va- If one takes the projected gradient of 
this expression one gets, after some manipulations, the 
relation 

Xa + {p + P)ZaEQ {Xa + Y^) + (W.M*' - h^)VbP = 0. 

( 3 ) 

We then rewrite Za as 

- h^)Xia. (4) 

This suggests introducing the covector 

Ca = = Daa - -DaP, (5) 

Hp + p) p 

which represents our non-perturbative generalization of 
the curvature perturbation on uniform density hypersur¬ 
faces or, more exactly, of its “spatial” gradient. Interest¬ 
ingly, the same linear combination has been introduced in 
the context of the long wavelength approximation, within 
a specific coordinate system, in Here, fa is defined 
geometrically and in the most general context. 

The spatial projection of Eq. 0 , via hj’, then yields 

e-“C(e“Cb)- = - C" + CVba) , ( 6 ) 

where we have used the familiar decomposition (see e.g. 
[ 3 ) VfeMa = aab+oJab + ^Qhab-UaUb with the (symmet¬ 
ric) shear tensor (Tab, and the (antisymmetric) vorticity 
tensor ujab', on the right hand side, the quantity 

^Ya - C^Xa = DaP -rT’aP (7) 

P 

represents the non-perturbative generalization of the 
nonadiabatic pressure and vanishes for purely adiabatic 
perturbations, i.e., when the pressure P is solely a func¬ 
tion of the density p. Note also that, for both Ca and 
Ta, one can replace the projected derivatives by ordinary 
partial derivatives respectively in Eqs. m and 0 . 

Equation is one of the main results of this work. 
It gives the fully non-linear evolution equation for the 


variable fa and corresponds to the non-perturbative gen¬ 
eralization of the familiar linear conservation law for adi¬ 
abatic perturbations (Ta = 0 ) on large scales, similar to 
the results recently obtained in the long wavelength limit 

i 

Introducing coordinates - To relate our approach to the 
more familiar coordinate-based perturbative approach, 
we now introduce a coordinate system, in which the per¬ 
turbed metric can be written as ds^ = af[—{l-\-2A)drf-\- 
2Bidx’'dr] {-jij S^ij) dx^dx^] (see 0). 

In the following, we decompose any quantity x in Ih® 
form x{v,P) = xijl) + + 51 ^X 2 (^. 2 ;*), where 

5xi and < 5 x 2 represent, respectively, the first and second 
order perturbations. Let us write down explicitly the 
components of our covector fa in this coordinate system. 
The zeroth order vanishes and, at first order, the spatial 
components are simply 

C«=5,Ci, fi=5a,-j5p,, ( 8 ) 

where a prime denotes the partial derivative with respect 
to the time coordinate 77 . To compute the component 
Co, it is useful to note that for any function x, one can 
write Dqx = uou^diX — u'^UidoXi where we have used the 
normalization of Since m* is first order, this implies 
that = 0 and • 

Expanding fi = ditt — {a/p)dip at second order, one 
finds, after some simple manipulations, 

cf ^ = 5,C2 + |r<5pi9.Ci', (9) 

C 2 = da2-^6p2-%5a\5p^+2^5p^5p\ + ]^ 5p\. 

P P P^ P \P J 

(10) 

Note that the conservation of fj'^^ and fj^'^ is equivalent 
to the conservation of (Ci and C 2 - 

We now relate a to the metric perturbations, consider¬ 
ing for simplicity only large scales. Neglecting gradients 
as well as first order vector and tensor perturbations, 
and writing if = —j'^^Sxij/S, one can use the equality 
3<i = Xau’^ to derive, up to second order, the relation 

a ~ In a— '(/> —'i/)^. ( 11 ) 

Thus, from Eq. (Cnj, we can relate C 2 to the quantity 
defined by Malik and Wands in |^, 

C 2 — C2(MW) - 2Ci(mW)' 

The expansion of the nonadiabatic term Tq = da P — 
{P/p)daP, which can be read from the expansion of fa 
by substituting P to a, enables us to write explicitly the 
conservation equation 0 at first and second orders. If 
one ignores the second term on the right hand side of 
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Eq. ®, one finds, at first and second order, respectively, 


c; - - 

C 2 ^ - 


n 


p + p 

n 

p + p 


ri, 


rEo- 


p 


EiCi', 


where the ri _2 are defined as the (’ 1 ^ 2 , and Ti. = a'ja. 
One thus recovers the results of very easily. Moreover, 
Eq. 0 generalizes them to any order in the perturbation 
theory and to all scales. 

Comoving “curvature” perturbation - In addition to 
the uniform energy density curvature perturbation, gen¬ 
eralized by our (^a, another useful quantity in the linear 
theory is the curvature perturbation on comoving hy¬ 
persurfaces, usually denoted by TZ. In our formulation, 
since the perturbed integrated expansion a replaces the 
usual spatial curvature and since our spatial gradients 
are defined with respect to the comoving observers, the 
generalization of TZ is simply 


are negligible, one can recover the familiar conservation 
law for adiabatic perturbations by noting that Xa^ which 
generalizes the comoving energy density perturbation of 
the linear theory can be shown to be negligible on 
large scales. Since Ca + TZa = — one then finds 

Ca — —TZa non-perturbatively. 

Note that TZa is not the only covariant quantity that 
can be connected, within a coordinate system, to the 
familiar comoving curvature perturbation of the linear 
theory. This is also the case for the projected gradient of 
the spatial scalar curvature (denoted by Ca in, e.g., 0 ), 
although there does not seem to be a simple relation at 
the non-perturbative level between TZa and Ca- 

Scalar field - It is also straightforward to extend our 
non-perturbative approach to a scalar field. Let us con¬ 
sider a scalar field (p with a potential Its energy- 

momentum tensor. Tab = dafidbfi — \gab {dcfid^cp + 2V), 
can be rewritten in the perfect fluid form by defining the 
four-velocity as the unit vector orthogonal to hyper¬ 
surfaces of constant p 0, 


TZa = -Daa = -Wa- (12) 

In the coordinate system introduced earlier, the spatial 
components are 


TZi = -dia - aa{vi -I- Bi), 


with u* = au*. At first order, 5a ~ —ip, and one rec¬ 
ognizes the familiar comoving curvature perturbation of 
the linear theory. 

It is relatively easy to derive the evolution equation 
for TZa- Equation can be seen as an equation for 
TZa = —Wa- Projecting orthogonally to u“, we find 

h^iZb+^&TZa = -7^'’(crba+W{,a)-^a^Mbd-iZa. (13) 

The term involving Za on the right hand side can be 
rewritten by using the shear constraint equation [see 
Eq. (4.17) of 0], 

-Za — + 0!%) + {aa + (^a)ub = Qa = 0, (14) 

where qa is the energy flux in the energy-momentum ten¬ 
sor, which, in our case, vanishes. Moreover, one can re¬ 
place iib using the Euler equation, iia = —Yal{p -I- P)- 
On rewriting Ya as Ta = Pa + c^Xa, and combining all 
this into Eq. CSJ, we finally obtain 


e-“C(e“7^b)- = 






3(p + p) 2ip + P)\ 


{Tb + clXb) 


-TZ^{aba + ^ba) - 




This is the fully non-linear evolution equation for TZa- 
As one can see, it is slightly more complicated than the 
evolution equation for (pa- However, when aab and ^ab 


Ua = -da(pl{-dc(pd‘'(py^^ = -dafilfi, (16) 

where the second equality (with = u“da(p) is a conse¬ 
quence of the definition of and the energy density and 
pressure as p = ^^/2 -|- V and P = — V- 

For a scalar field, the comoving curvature perturba¬ 
tion is particularly convenient because the (p = const 
hypersurfaces are orthogonal to defined in Since 
Da(p = 0, one sees that TZa = —Daa can be reexpressed 
as 


TZa = -daa + ^da(p, 

and one recognizes the comoving curvature perturbation 
of the linear theory. The generalization of the Mukhanov 
variable Q for the quantization of the scalar field-gravity 
system in the non-linear theory is thus Va = e°‘{(pl dpTZa- 
In a coordinate system, TZa can be easily expanded to 
first and second orders to make contact with previous 
results. As for (pa^ one finds, at first and second orders, 
expressions of the form 

= 9 , 7^2 + ^ScPid,TZ[, 

9 

where the quantities TZi and TZ2 can be deduced from 
Cl and C 2 , given in 0 and 0, by changing the overall 
sign and replacing p by (p- When first order vector and 
tensor modes can be neglected, up to small scale terms 
and an overall sigm TZ = TZi + DjZb coincides with the 
variables used in (see also [ig ) either in the uniform 
field gauge, which can be obtained by setting 5(p = 0, or 
the uniform curvature gauge, obtained by setting 5a = 0. 
On writing a in terms of the metric components, and 
taking the large scale limit, one can also show that TZ 
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can be simply related to the conserved quantities defined 

in 110 . 

For a scalar field {cOab = 0)) if <^ab can be neglected, one 
can also generalize the well known first order property 
that the nonadiabatic term Fq is negligible in the large 
scale limit. This result, which extends what is known up 
to second order 0 , is obtained by noting that DaV = 0 
(since Dacj) = 0), which implies F^ = (1 — Cg)Xa ~ 0. 
In this case, for a single scalar field, the right hand 
side of Eq. m is negligible and TZa is conserved, on 
super-Hubble scales, to all orders in perturbation the¬ 
ory. If more fields are present, then F^ cs —2DaV 0 
and the non-linear perturbation TZa, and thus the non- 
Gaussianity, can evolve on large scales during inflation. 

Conclusion - In this Letter, we have defined, respec¬ 
tively in © and non-perturbative generalizations 

of the so called uniform energy density curvature pertur¬ 
bation and comoving curvature perturbation, using the 
integrated expansion a, which represents the number of 
e-folds of the local expansion measured by an observer 
comoving with the fluid. 

Although our perturbed quantities are similar to those 
advocated in recent works based on the long wavelength 
approximation their originality is that they are 

constructed geometrically, independently of any coordi¬ 
nate system. Furthermore, we have been able to derive 
evolution equations governing these quantities, which are 
fully non-linear and exact at all scales. In particular, 
Eq. obtained directly from the local conservation 
of the energy-momentum tensor and independently of 
Einstein’s equations, is valid in any metric theory of 
gravity, including scalar-tensor theories or induced four¬ 
dimensional gravity in brane-world scenarios. Thus, our 
relations generalize the familiar conservation laws of the 
linear theory for adiabatic perturbations. 

Another advantage of our formulation is that it is 
rather easy to relate to the coordinate-based perturbative 
formulation. Interestingly, our extremely simple non¬ 
linear expressions, like automatically contain the 
seemingly complicated second-order structure that has 
been obtained in the coordinate-based approach, which 
we can easily rederive and extend to higher orders. 

Finally, let us stress that our formulation, since it keeps 
track of all small scale terms, may simplify the study of 
non-linear perturbations even inside the Hubble radius. 
This would be very useful, in particular to compute the 


non-Gaussianity from inflation [13 or to study the back- 
reaction of small scale modes on larger scales. 

We thus believe that this approach might represent a 
significant contribution to the study of primordial cos¬ 
mological perturbations beyond the linear order. 
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